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Abstract 

This note provides a general construction, and gives a concrete example of, forced ordinary 
differential equation systems that have these two properties: (a) for each constant input u, all 
solutions converge to a steady state but (b) for some periodic inputs, the system has arbitrary (for 
example, "chaotic") behavior. An alternative example has the property that all solutions converge 
to the same state (independently of initial conditions as well as input, so long as it is constant). 

1 The question 

We consider systems 

x{t)^ f{x{t),u{t)) 

of forced ordinary differential equations, where u — u{t) is an input, assumed scalar- valued for simplicity. 

It is well-known fact that a system that exhibits periodic trajectories when unforced may well exhibit 
chaotic behavior when forced by a periodic signal, as shown by Van der Pol equations with nonlinear 
damping, and nonlinear-force Dufhng systems. However, the following fact seems to be less appreciated. 

Assume that the system has the property that, whenever u{t) is constant (any value), all solutions 
converge to a steady state (which might depend on the input and on the initial conditions). Is it possible 
for such a system to behave pathologically with respect to periodic inputs (no "entrainment")? We 
outline a general construction together with a specific example of this phenomenon. 

A second construction is also given, which has the property that all solutions for constant inputs converge 
to the same state (independently of initial conditions as well as the input value, so long as it is constant). 
This second construction displays a more interesting phenomenon, but it is only done for an example 
(Lorentz attractor), and it seems difficult to prove rigorously that its solutions when forced by a periodic 
input are truly chaotic. 



2 The general construction 

The general construction is illustrated in Figure [TJ 

The block called "filter 1" consists of a linear system x = Ax + Bu with output y = Cx + Du having 
the property that its transfer function W{s) (the Laplace transform of the impulse response Ce^^B) is 
not identically zero but has a zero at the origin: W{0) = 0. 
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Figure 1: General construction 



The difFerentiable function a : M ^ M>o has the property that a(0) = but w{t) — a{y{t)) saturates at 
the value 1 if y{t) is outside a small neighborhood of zero. 

The block called "filter 2" consists of a scalar linear system p = —p + w. If the signal w{t) converges to 
zero, then p{t) also converges to zero. However, if the signal w{t) « 1 for a large set of times t, then p{t) 
will, after a transient, stay always positive and bounded away from zero, through an averaging effect. 

Finally, the block named "system" has the form z(t) = p{t)f{z{t)), where / is an arbitrary vector field 
with bounded solutions. Provided that the scalar signal p(t) remains positive and bounded away from 
zero, the solutions of the system will be time-rescaled versions of solutions of i = f{z), and hence have 
the same pathologies. 

When the input u{t) is constant, the fact that the transfer function vanishes at zero means that y{t) — > 
as i — > oo. This implies that w{t) 0, and so p{t) ^ as well. Thus the velocity z{t) 0, and 
solutions become asymptotically constant. 

On the other hand, suppose that the input is u{t) = sinuit, where W{iu!) ^ 0. This means that the 
output will be asymptotically of the form ylsin(wt + (p), for some phase £ R and some amplitude 
A = \W{iuj)\ ^ 0. Assuming that A is large enough, passing y{t) through a will result in a signal w 
that is mostly saturated near the value 1, so p(t) will stay away from zero and the solutions of the z 
subsystem will look (asymptotically) like those of i = f{z). 



3 An explicit example 

Let us take "filter 1" to be the system x = —x — u with output y = x + u, the function a{y) = y"^ / {K + y'^) 
with K = 1/10, and the vector field corresponding to the Lorentz attractor with parameters in the chaotic 
regime. Specifically, the complete system is as follows. 

X = —x — u 

p = —p + a{x + u) 

i = P{s{i^~0) 

i{j = p{r^ _ ^ _ 

C = pi^i^-bC) 

with Prandtl number s = 10, Rayleigh number r — 28, and b ~ 8/3. Figure [2] shows typical solutions of 
this system with a periodic and constant input respectively. 



4 An alternative example 

Instead of multiplying the Lorentz vector field by p(t), we may use p{t) to eliminate the multiplicative 
terms, making the system asymptotically linear and time-invariant, and also changing one of the coeffi- 
cients in such a way that the corresponding linear system is globally asymptotically stable. Simulations 
(using a Montecarlo sampling of initial conditions as well as input magnitudes) indicate that the system 
remains chaotic when forced by sini, but all solutions converge to zero when u{t) is constant. Specifically, 
the complete system is as follows. 

X = —X — u 
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Figure 2: Some solutions of example system. Initial conditions are x{0) = 5, ^(0) = 1, and all remaining 
variables zero. Green: inputs are u{t) = s'mt (left panel) and u{t) = 10 (right panel). Blue: £^{t). Note 
chaotic- like behavior in response to periodic input, but steady state in response to constant input. 



p = —p + a{x + u) 

i = 10(^-0 

C = pe^-(8/3)C 

Observe that when p = Q one has, for the last three variables, the stable linear system: 



io(V-0 



and when p — 1 the system is the usual Lorentz attractor 

^ - io(V^-0 

C = ?^-(8/3)C. 

In order to achieve p « 1 for the input sint, we now pick K = 0.0001 in the definition of a. Figure [3] 
shows typical solutions of this system with a periodic and constant input respectively. The function 
"rand" was used in MATLAB to produce random values in the range [—10, 10]. To better understand 
the behavior of this system. Figure [3] shows the value of the input and the variable p{t), for the same 
input and initial conditions as in Figure [31 



5 Remark 

One could clearly extend the construction to "interpolate" between any two dynamics: using z{t) = 
p{t)fi{z{t)) + (1 — p{t))fo{z{t)) means that, for constant inputs, asymptotically the system wiU behave 
like i = fo{z{t)), while for u{t) = sint it will behave like i = fi{z{t)). Establishing a precise result 
in this generality would require a detailed study of chain-recurrent sets and other objects associated to 
asymptotically autonomous dynamics. We leave that to future work. 
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Figure 3: Simulation of second example, done with the following randomly-chosen input and initial 
conditions: = 1.89, a;(0) = 2.95 p(0) = -0.98, ^(0) = 0.94, V(0) = -4.07, C(0) = 4.89. Green: 
inputs are u{t) = sint (left panel) and u{t) = 5.13 (randomly picked, right panel). Blue: £^{t). Note 
chaotic-like behavior in response to periodic input, but steady state in response to constant input. 
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Figure 4: Inputs and variable p{t) for simulation shown in Figure 31 Note that p{t) « 1 and p{t) « for 
i>0. 
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